Clustering of heavy particles in the inertial range of turbulence 
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A statistical description of heavy particles suspended in incompressible rough self-similar flows is 
developed. It is shown that, differently from smooth flows, particles do not form fractal clusters. 
They rather distribute inhomogeneously with a statistics that only depends on a local Stokes number, 
given by the ratio between the particles' response time and the turnover time associated to the 
observation scale. Particle clustering is reduced when increasing the fluid roughness. Heuristic 
arguments supported by numerics are used to explain this effect in terms of the algebraic tails of 
the probability density function of the velocity difference between two particles. These tails are a 
signature of events during which particle couples approach each other very closely. 

PACS numbers: 47.27.-i, 47.51. +a, 47.55.-t 
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In turbulent flow, finite-size particles heavier than the 
carrier fluid, as e. g. water droplets in air, possess inertia 
and have a finite response time to the fluid motion. Thus 
their dynamics markedly differs from that of simple trac- 
ers and, in particular, such inertial particles distribute in 
a strongly inhomogeneous manner even if the underlying 
flow is incompressible. Modelling these fluctuations in 
the particle concentration is important in engineering [l| , 
cloud physics Q , and planetology . 

The turbulent motion of the carrier fluid spans many 
active spatial and temporal scales 0. Below the Kol- 
mogorov length-scale, viscous dissipation dominates; the 
velocity field is diffcrcntiable and is characterized by a 
single time scale. There the motion of inertial particles 
is governed by the fluid strain and their dissipative dy- 
namics leads their trajectories to converge to a dynam- 
ically evolving attractor. For any given response time 
of the particles, their mass distribution is singular and 
generically scale-invariant with multifractal properties at 
small scales 0, IS 0- Above the Kolmogorov scale, the 
fluid velocity field is not smooth anymore but, accord- 
ing to the Kolmogorov 1941 theory, almost self-similar 
with Holder exponent h = 1/3 Q . This so-called inertial 
range is characterized by a broad-spectrum of time scales. 
However, the finite response time of the suspended par- 
ticles introduces a new scale, breaking self-similarity in 
the particle distribution. This is consistent with the ob- 
servation that particles typically concentrate on different 
scales with the largest deviation from uniformity aris- 
ing when their response time is of the order of the eddy 
turnover time [1, EJ. As already noticed in Q, these 
deviations are not scale-invariant and have a different 
origin from those observed in the viscous range of turbu- 
lence. With few exceptions |ToL ITlL fl~3 | . clustering in the 
inertial range received considerably less attention than 
small-scale clustering. 

In this Letter we focus on second-order statistics of the 



particle distribution at scales within the inertial range. 
These statistics can be completely described in terms of 
the relative motion of particle pairs. Within the inertial 
range, two concurrent mechanisms responsible for par- 
ticle clustering can be identified: a dissipative dynam- 
ics due to their viscous drag with the fluid and ejection 
from persistent vortical regions by centrifugal forces . 
By modelling the carrier flow as a rough, self-similar 
random velocity field which is (5-correlated in time, we 
eliminate the latter effect: the absence of any persistent 
structure in the considered carrier flow ensures that cen- 
trifugal forces play no role. This model pertains to very 
heavy particles whose response time is much lar ger than 
the typical correlation time of the ambient fluid [lj, njj . 
We show heuristically and numerically that the scale- 
invariance of the velocity field does not extend to the 
particle distribution, and that roughness of the carrier 
velocity weakens clustering. This effect is explained by 
the dependence of the relative velocity distribution on 
the fluid velocity Holder exponent. 

The relative motion of two particles is described by 
their separation R that obeys the equation 0, ^| 

R = — \R-5u(R,t)] , (1) 
r L J 

where the dots denote time derivatives, r is the Stokes 
time, and Su(r, t) = u(x + r,t) — u(x, t) is the fluid ve- 
locity difference. The velocity u is assumed to be a sta- 
tionary, homogeneous and isotropic Gaussian field with 
correlation 

( Ui (x, t) Uj (x', t')) = 2D Q S tJ - B^x - x') S(t - 1') , (2) 

where D is the velocity variance. For rough self-similar 
flows, the function B takes the form 

Bij (r) = D 1 r 2h [{d - 1 + 2h) 8 VJ - 2h r.^/r 2 ] , (3) 

where r = \r\, d is the space dimension, h 6 [0,1] the 
Holder exponent of the carrier velocity field, and D\ a 
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S(r) S(r) 
FIG. 1: (a) Local correlation dimension S(r) versus the scale-dependent Stokes number S(r) = Dir/r 2 ' 1- ' 1 ' for various values 
of the Holder exponent h in two dimensions d = 2. (b) Ratio between the local exponent 7(r) of the particle velocity and h 
versus S(r). Different symbols (colors online) refer to different values of the particle response time t. 



constant measuring the intensity of turbulence. This 
kind of velocity field was introduced by Kraichnan [lf3 | 
to model passive scalar transport. 

By defining s — t/r and rcscaling R by the obser- 
vation scale £, it is easily seen that the above dynam- 
ics, and thus all the statistical properties of particle 
pairs at scale £ only depend on the local Stokes num- 
ber S(£) = D\T I '^ 2 ( 1 -' 1 ). This dimensionless quantity, 
first introduced in 0, is the ratio between the particle 
response time r and the turnover time at scale t. It 
measures the scale-dependent effects of inertia. At large 
scales (£ — > oo) inertia becomes negligible (S(£) — > 0) and 
particles recover the incompressible dynamics of tracers. 
Conversely since S(£) — * oo for I — ► 0, inertia effects 
become dominant at small scales and the dynamics ap- 
proaches that of free particles. For both very large and 
very small values of the Stokes number, the particles dis- 
tribute uniformly in space. Strong inhomogeneities ap- 
pear for S(£) ~ 1. 

Note that in unbounded carrier flows, the separation 
between two particles asymptotically grow indefinitely 
with time and thus the dynamics Q never reaches a 
statistical steady state. However, real turbulent flows 
are bounded, allowing for statistical equilibrium to be 
reached. Boundary conditions are thus implemented in 
the considered model by imposing, for instance, reflec- 
tion of the inter-particle distance at \R\ = L. Clearly, 
since self-similarity is broken by the presence of bound- 
aries, the aforementioned scaling arguments apply only 
at scales £ <C L. 

For smooth carrier flows (h = 1), there is a unique time 
scale so that the dynamics only depends on the global 
Stokes number S = D\t. Inhomogeneities in the particle 
distribution can be quantified by d — 2?2, where T>2 is the 
correlation dimension 

V 2 = lim 8{r), S(r) = d (lnP 2 (r)) /d (lnr) , (4) 



P2{r) denoting the probability that \R\ < r. In 5- 
correlated smooth flows, just as in real suspensions, the 
correlation dimension T>2 non-trivially depends on the 
Stokes number [l^j . 

For non-smooth but Holder-continuous flows, T>2 = d 
for all particle response times. Information on the in- 
homogeneities of the particle distribution is entailed in 
the local correlation dimension S(r) defined in (@J). From 
above arguments, S(r) is expected to depend only on h 
and on the local Stokes number S(r) when r <C L. This 
is confirmed numerically for d = 2 as shown in Fig. ^i, 
where the local dimension S(r) is represented as a func- 
tion of S(r) . The collapse of the data for various response 
times r demonstrates the dependence on the local Stokes 
number. Comparing various values of the exponent h, we 
observe that when the fluid becomes rougher, the inten- 
sity of clustering weakens. In particular, the minimum of 
8(r) gets closer to d = 2 as h decreases. Notice that for 
h = 1, the Stokes number S does not depend on r and 
data refer to the actual value of the correlation dimen- 
sion, which is well defined for smooth flows (see for 
details). 

We now turn to the typical velocity difference R be- 
tween two particles and its dependency on the separation 
R. For smooth flows, when \R\ — > an algebraic behav- 
ior of the form \R\ ~ |i?| 7 is observed, defining a Holder 
exponent 7 for the particle velocities. This exponent de- 
creases from 7 = h = 1 for S = 0, corresponding to a 
differentiable particle velocity field, to 7 = for S — > 00, 
that means particle moving with uncorrelated velocities 
|l4j . In non-smooth flows the exponent 7 is asymptoti- 
cally equal to the fluid Holder exponent h at large scales 
(S(r) — > 0), while it approaches at very small scales 
(iS(r) — ► 00). Therefore, similarly to S(r), all relevant in- 
formation is entailed in the scale dependence of the local 
exponent "f(r). As for the local correlation dimension, 
this exponent only depends on the fluid Holder expo- 
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nent and on the local Stokes number; this is confirmed in 
Fig. ^d, showing the ratio j(r)/h versus S(r) for various 
values of h. It is worth noticing that the transition from 
7(7") = h to j(r) = shifts towards larger values of the 
local Stokes number and broadens as h decreases. The 
fact that j(r) = h for r — > 00 implies that the particles 
should asymptotically experience Richardson diffusion 0| 
just as simple tracers. 

To get a deeper understanding of the mechanisms of 
clustering, we transform the equations of motion Q)-© 
into a system of stochastic differential equations with ad- 
ditive noise. Adapting to rough flows the strategy first 
proposed in |l7l | for smooth carrier velocities fields and 
used in [J^, we make the following change of variables: 



X = {t/L 2 )(\R\/L)-( 1+ V R-R, 
Y = {T/L 2 ){\R\/L)- (1+h) \RAR\ 
Z = (| J R|/L) 1 -\ 



(5) 



The variables X and Y refer to the longitudinal and 
the transverse dimensionless velocity differences, respec- 
tively. When d = 2, the dynamics Q in the four- 
dimensional phase-space (R, R) reduces to 



X = -X- 
Y = -Y-(l + h)Z 
Z = (1 - h) X . 



Z~ l (hX 2 ■ 

-1 



-Y 2 ) 
XY 



f 771 (s) , 



(G) 
(7) 
(8) 



Now the dots denote derivatives with respect to s = 
t/r; rj\ and 772 are independent white noises with vari- 
ances 2S{L) and 2 (1 + 2h)S(L), respectively; S(L) = 
D\T / L 2 ^ h ^ is the Stokes number associated with the 
system size. Periodic boundary conditions in physical 
space amount to considering reflective boundary condi- 
tions at Z = 1 ; Y is ensured to remain positive by reflec- 
tive boundary conditions at Y = 0. Rcscaling \R\ with 
A, and thus Z with A 1 ~ h leads to transform X and Y to 
\ 1 ~ h X and A 1 ~ h Y in order to confine the scaling factor 
in the noise. This again amounts to considering the same 
dynamics with a scale-dependent Stokes number S(XL). 
The system (JBJ-lJSJ can be efficiently implemented nu- 
merically and was used to produce the data described in 
this Letter. 

Figure [2 sketches the dynamics in the (X, Y, Z) space. 
The line X = Y = (its physical meaning is that the 
particles have relaxed to the same velocity staying at an 
arbitrary distance) acts as a stable fixed line for the drift 
terms in A typical trajectory spends a long time 

diffusing around this line, until the noise realization be- 
comes strong enough to escape from its neighborhood. 
When this happens with X > 0, the quadratic terms in 
the drift drive the trajectory back to the stable line. On 
the contrary, if X < and hX 2 +XZ-Y 2 < 0, the drift ac- 
celerates the trajectory towards larger negative values of 
X. Then Z decreases until the quadratic terms in ©-0 
become dominant. The trajectory then loops back in the 
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FIG. 2: Phase-space picture of the system (JSJ-JSJ for h — 
7/10. The thin smooth lines represent effects of the drift. A 
random trajectory of the system with S(L) = 1 is shown as a 
bold line; it performs a large loop from X < to X > 0. (a) 
The full (X, Y, Z) space, (b) and (c) projections in the Z = 
and Y — planes, respectively. 

{X , y)-plane, approaching the stable line from its right. 
These loops are the events during which Z (and hence the 
inter-particle distance R) becomes substantially small. 

As we now show, these loops are responsible for power- 
law tails in the probability density function (PDF) of the 
dimensionless velocity differences X and Y. This behav- 
ior can be understood by extending to the rough case the 
arguments developed in 0] for smooth flows. For this we 
consider the cumulative probability P < (x) = Pr {X < x) 
with x <C — 1 , which can be estimated as the product of 
(i) the probability to start a sufficiently large loop that 
reaches values more negative than x and (it) the fraction 
of time spent by the trajectory at X < x. To estimate 
these two contributions, we assume that within a distance 
of order unity from the line X = Y = the quadratic 
terms in the drift are negligible, so that X and Y be- 
have as two independent Ornstein-Uhlcnbcck processes. 
At larger distances we retain only the quadratic terms 
responsible for the loops. 

Within this simplified dynamics, a loop is initiated 
at a time sq for which X$ — X(sq) < — 1 and Yq = 
y( s o) ^ l-^o I- The maximum distance from the sta- 
ble line is attained at a time s* for which X(s*) is of 
the order of -Y(s*), that is when X(s*) = -pY(s*), 
(3 being an arbitrary constant. When the noise is ne- 
glected, one straightforwardly obtains: Y(s)/X(s) = 
Y Z /{X a Z + (l-e s °- s ) [A 2 + y 2 ]} and the radius of 
the loop can be estimated as 



LY( S *)| 



[X + Z Q 



yh 



Yn 



[ 1+/3 2](l + h)/2 



(9) 



For the trajectory to reach values X < x <C — 1, the ra- 
dius has to be larger than —x, and thus Y(so) has to be 
smaller than |x| _1/ '' 1 . The joint PDF of X and Y at time 
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FIG. 3: PDF of X in log-log coordinates for S(L) = 1 and 
various values of h. In all cases, power-law tails are observed. 
Inset: exponent a of the algebraic tail as a function of the 
fluid velocity Holder exponent h; the theoretical prediction is 
represented as a dotted line. 

So is given by the dynamics close to the line X = Y = 0. 
As it is finite for Y — >0, the contribution (i) can be es- 
timated to be oc The contribution (ii) can be 
obtained as follows. Far from the stable line, the dynam- 
ics can be approximated by the deterministic part, hence 
the fraction of time spent at X < x is oc Yq oc |ar| -1 /'*. 
Putting together the two contributions, we obtain that 
P < (x) oc |x| -2 /' 1 when x <C — 1. The negative tail of the 
PDF of the longitudinal velocity difference X behaves as 
a power-law oc \x\~ a with a = 1 + 2/h. This gives a = 3 
for smooth flows (h = 1) as previously derived [Til ]. Dur- 
ing the large loops, the trajectories equally reach large 
positive values of X and of Y. This gives again a frac- 
tion of time x~ x l h spent at both X and Y larger than 
x ^> 1. Hence, the PDF of both the longitudinal and 
the transversal velocity differences have algebraic left and 
right tails. 

As shown in Fig. |3 the presence of power-law tails in 
the PDF is confirmed numerically, with perfect agree- 
ment between the measured values of a and the predic- 
tion a = 1 + 2/h (see inset). Let us comment on the h 
dependence of this exponent a. The large loops respon- 
sible for the algebraic tails correspond to events in which 
particles approach each other very closely; they are the 
basic mechanism of particle clustering. The probability 
to enter such loops decreases significantly when h — * 0. 
Moreover, it is straightforward to check from I©-© that 
during the loops Z(s) oc Zjj when Zq <C 1. Hence it gets 
less and less probable to reach smaller values of Z as 
h decreases. Combined together, these two effects ex- 
plain why particle clustering is weakened in rough veloc- 
ity fields and why it is more efficient in smooth flows. 

The change of variables JSJ) can be equally applied in 
three dimensions, leading to a dynamics different from 
©-©■ Therefore understanding to what extent the 



above findings extend to the 3D case remains an open 
question; work in this direction is under development. 

To conclude, let us comment on the implications of 
this work on the study of heavy particles in realistic tur- 
bulent flows. There, particle clustering is simultaneously 
due to ejection from eddies and to a dissipative dynam- 
ics. The considered model flow isolates the latter effect. 
It is probable that power-law tails for velocity differences 
can be present in realistic settings as well. However, it 
is not clear if the results on clustering are affected by 
the presence of persistent structures: particle ejection 
from eddies may form voids and thus v ery strong inho- 
mogeneities in the particle distribution |lfj |. This could 
overtake dissipative-dynamics mechanisms. 
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